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Abstract While very often a hot intergalactic medium (IGM) is optically thin to 
continuum radiation, the optical depth in resonant lines can be of order unity or larger. 
Resonant scattering in the brightest X-ray emission lines can cause distortions in the 
surface brightness distribution, spurious variations in the abundance of heavy elements, 
changes in line spectral shapes and even polarization of line emission. The magnitude 
of these effects not only depends on the density, temperature and ionization state of 
the gas, but is also sensitive to the characteristics of the gas velocity field. This opens a 
possibility to use resonant scattering as a convenient and powerful tool to study IGM 
properties. We discuss the application of these effects to galaxy clusters. 

Oh' Keywords X-rays • clusters of galaxies • radiative transfer 
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CZ) 1 Introduction 

Hot (10 — 10 K) X-ray emitting astrophysical plasmas are often optically thin with 

respect to free-free absorption. One famous exam ple is the coronae of the Sun and 
^ ' other stars. It was realized early (e.g. lElwerdll95a ) that in the solar corona the optical 

depth in certain resonant lines of oxygen, neon and iron can be larger than unity and a 
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radiative transfer problem has to be s olved in order to get the correct line ratios (e.g. 
lActonlll978l : iRueee fc McKenzielll985l ). 

Below we concentrate on the resonant scattering in the hot plasma of galaxy clus- 
ters. Hot plasma is the dominant constituent of baryons in clusters, making 10-15% of 
the total cluster mass of 10 —10 Mq. Bremsstrahlung, recombination continuum 
and emission lines excited by electron collisions are the main emission components in 
this plasma. In spite of such a large gas mass, the optical depth of the hot plasma 
for free-free absorption is extremely low in the X-ray regime. This is immediately 
clear from a comparison of the observed X-ray luminosity of the brightest clusters, 
Lx ~ 10 erg s _ , with the black body radiation of a Mpc size object with a tem- 
perature T ~ 10 7 - 10 8 K: 



- L * -i k in-33 ( L x \(R\ ( T y 4 
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where R is the characteristic cluster size. The depth for Thomson scattering is also 
low: 



^ = ^R = 2 ■10- 3 ( IF ^ FI ) (^), (2) 

where n e is the electron density. This means that effectively clusters are transparent 
for X-rays. As in the solar corona, this statement is n ot valid for the brightest X-ray 
resona nt lines of the most abundant elements like iron. iGilfanov. Sunvaev fc Churazovl 
(|1987l ) have shown that the optical depth in the 6.7 keV line of He-like iron in the 
brightest cluster A426 (Perseus) is larger than unitjU and therefore effects of reso- 
nant scattering should be visible in the surface brightness distribution, element abun- 
dance variations or in the spectral shape of the line. It was immediately clear that 
all these effects are sensitive to the turbulent broadening of the line, making reso- 
nant scattering a unique tool to probe the characteristic amplitude of the gas ve- 
locities using telescopes with limited energy resolution. The impact of resonant scat- 
tering on surface brightness profiles has been studied via radiative transfer calcula- 
tions (e.g. Gilfanov. Sunvaev fc Churazovlll987i : IShigevam a 1998; C hurazov et alj|2004 : 



I zhuravW et ab l ^Ohh and searched for ^s^v^o^f^ .rs. grouns of ga la^s 
and individual massive elliptical galaxies. 

It has been realized that the scattering should also lead to a pola rization of the 



it has been realized that the scattering should also lead to a polarization ol the 
same lines ( Sazonov. Churazov fc SunvaevN2002l ; IZhuravleva et alj|20ich . with the po- 



larization degree varying from ~ at the center of a cluster to tens of per cent in the 
cluster outskirts (although the unpolarized continuum and nearby unpolarized lines 
will decrease this figure for an instrument with modest spectral resolution). 

There have been numerous experimental studies aimed at detecting effects of reso- 
nant scattering in rich clusters and elliptical galaxies. For the Pers eus cluster, an initial 
analysis s uggested that resonant scattering effect s were present l|Molen di et al.||l998: 



Akimoto et alJl999f ). although subsequent s tudies (Churazov et al. 2004; Gastaldcll o fc Molendil 



2004J) did not suppo rt these results (see also lEzawa et al.ll200ll ; lDupke fc Arnaud 



2001 



Sanders et al.ll2004r ) . Several other clusters , including M87, Centaurus and A2199, were 
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searched for resonant scattering effects (c.g Kaastra, Blocker, & Mcwc 1999; Akimoto ct al. 



19991 .12000: Sand ers fc Fabianl20 06: Mat hews. Buote. fc Brig henti 2001; Sakelli ou et al 



For pure thermal broadening of the line. 



12002 ) . with mixed results. For M87 and Perseus, resonant scatt ering was mentioned as 



one of the possible explanations of the central abundance dip i Bohringer et al.ll200ll ; 
ISanders et al"]|2004) . 



In one objec t - the elliptical galaxy NGC4636 - m ost studies ( Xu et al.l 2002: 



iKahn et al.ll2003l ; IWerner et al.ll2009l ; lHavashi et al.l l2009) agree that resonant scatter- 
ing modifies the flux ratio of the Fe XVII lines at 15.01 A and 17.05 A. The same 
study also found evidence, albeit at low statistical confidence, for resonant scattering 
in several other systems: NGC1404, NGC5813 and NGC4472. 

In addition to the above mentioned studies several o ther ideas related t o reso - 
nant sca ttering in hot cluster p lasma have been proposed. Kroli k fe Raymond! ( 198a ') 



(see also Sarazin 1989; Mo lnar. Birkinshaw. fe Mushotzkvll2006r ) suggested to use mea- 
surements of the optical depth of the cluster plasma in lines together with the X-ray 
surface brightness to determine both the gas density and the linear size of the cluster. 
The known linear size can then be used for a cosmological distance-angular diameter 
test. lSazonov. Sunvaev. fe Cramphornl (2002J) pointed out that scattering of continuum 
emission in a resonant line can be used to constrain the past X-ray luminosity of AGN 
in clusters of galaxies. Effects of resonant scattering of the cosmic X-ray background 
(CXB) were sugg ested to play a role in the line emissivity of the WHIM - warm-hot 
intergalactic gas i Churazov et al.l I2001T I. provided that one can resolve a substantial 



fraction of point sources contributing to the CXB. 

The structure of this paper is as follows: in section I2TT1 we discuss the most promis- 
ing resonant lines and calculate the optical depths for several objects; the impact of 
resonant scattering on line surface brightness profiles and on abundance measurements 
is discussed in section [3] in section [4] we analyze the sensitivity of the resonant scat- 
tering effects to the gas velocity field; in section [5] the polarization of X-ray lines is 
discussed; in section|S]we illustrate the effects of resonant scattering on the line shape; 
in section [7] we briefly discuss the possible uses of resonant scattering for cosmological 
tests f ^7.1|l . for searches of powerful AGN outbursts in the past (" 37. 2[) and for studying 
the WHIM (373J- 



2 Basics of resonant scattering 

2.1 Optical depth 

The cross section for scattering at the center of a resonant line can be written as 

° = ^E^> (3) 

where r e is the classical electron radius and / is the absorption oscillator strength of 
a given atomic transition and AE]j is the Doppler width of the line. In a plasma with 
a temperature typical of galaxy clusters, the line width is determined by the velocities 
of thermal and turbulent motions, rather than by the radiative width. For example, 
for the 6.7 keV line of He-like iron the radiative width is ~ 0.3 eV, while thermal 
broadening is ~ 3 eV for a 5 keV gas. The Doppler width of the line is defined as 

AE D = E JM^ + ^fH , (4) 



where A is the atomic mass of the corresponding element, m p is the proton mass and 
Vturb is the characteristic turbulent velocity. Vturb is often parametrized as Vt ur b = 
c s M, where M is the Mach numbeiQ and the sound speed in the plasma is c s = 
■\/~fkT/fimp, where 7 = 5/3 is the adiabatic index for an ideal mono-atomic gas and 
jj, = 0.61 is the particle mean atomic weight. We can rewrite the previous expression 
as 



AE D = E 
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The optical depth of the cluster at the center of the line is then r — I agriidl, where I 



I 



is the distance along the photon propagation direction and n; is the number density of 
ions in the ground state of a given transition. For cluster conditions all ions are in the 
ground state (the frequency of collisions or any other excitation process are negligible 
compared to the radiative decay rate of the excited state) and n, = n p Z8i, where n p 
is the density of protons, Z the abundance of the element relative to hydrogen and <5j 
is the fraction of the element in the appropriate ionization state. That is 
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(6) 



Obviously, to have a large optical depth in the line one needs an astrophysically abun- 
dant element, a large oscillator strength of the transition, an appreciable ionization 
fraction of a given ion and a small line width. If the line width is dominated by ther- 
mal broadening then the lines of the heaviest elements have an advantage over lighter 
ele ments, since the thermal broadening sc ales as 1/VA. 

iGilfanov. Sunvaev fc Churazovl (|1987l ) give a convenient expression for the optical 
depth (from the cluster center to the observer) at the midpoint of a Doppler-broadened 
resonance line when the isothermal gas density distribution can be approximated by a 
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where r is the gamma- function; parameters of the 6.7 keV line of FeXXV are used; and 
Z/Zq is the abundance of a given element relative to the solar abundance of iron. In 
Table[TJwe quote the optical depths of several promising lines for 3 objects, representing 
different mass limits - from an elliptical galaxy (NGC4636) to a rich cluster (Perseus). 
The optical depths were calculated using the observed temperature and density radial 
profiles. Pure thermal broadening of the lines was assumed. 

From Table [T] it is clear that different lines can be optically thick in objects having 
drastically different masses and temperatures. For objects with T > 3 keV (e.g. Perseus 



2 This form of parametrization is rather arbitrary. For isotropic turbulence in a gas with 
7 = 5/3, the value of M = 1 corresponds to the energy in turbulent motions being equal to 



0.83 times the thermal energy density 



1 



7 - 1 /BTip 
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Table 1 Oscillator strength, Rayleigh scattering weights and optical depth of the most promis- 
ing X-ray lines for the elliptical galaxy NGC 4636 and the M87/Virgo and Perseus (A426) 
clusters. The optical depths were calculated using the observed temperature and density pro- 
files and assuming the fla t abundance profile Z = 0.5Zq (relative to the Solar abundances of 
Anders & Grcvcssc (1989)), collisional ionization equilibrium and pure thermal broadening of 



the lines. 














Ion 


E, keV 


/ 


U> 2 


r, NGC 4636 


t, Virgo/M87 


r, Perseus 


O VIII 


0.65 


0.28 


0.5 


1.2 


0.34 


0.19 


Fe XVIII 


0.87 


0.57 


0.32 


1.3 


0.0007 


1.510" 7 


Fe XVII 


0.83 


2.73 


1 


8.8 


0.0005 


2.8-10" 8 


Fc XXIII 


1.129 


0.43 


1 


0.016 


1.03 


0.16 


Fe XXIV 


1.168 


0.245 


0.5 


0.002 


1.12 


0.73 


Fe XXV 


6.7 


0.78 


1 


0.0002 


1.44 


2.77 



cluster) the He-like iron line at 6.7 keV is a likely winner; while for cool systems (e.g. 
NGC4636 with T ~ 0.6 keV) the Ne-like iron line at 0.83 keV has the largest optical 
depth. 



2.2 Phase function 



For any transition, the resonant scattering can be represented as a combination of two 
processes: iso tropic scattering with a weight u>i an d Rayleigh scattering with a weight 
u>2 — 1 — Wi I Hamiltonlll947l : IChandrasekharjll950r ) . The weights w\ and u>2 depend on 
the total angular momentum j of the ground level and on the difference between the 
total angular momenta of th e excited and gro und levels Aj (=±1 or 0). The expressions 
for the weights are given bv lHamiltonl 1 1947T ). For a subset of the most promising lines 
the W2 values are given in Table [T] 

If the radiation is initially unpolarized, then the probability for a photon to scatter 
into a unit sold angle at an angle 6 with respect to the incident direction is 



Pfa) 
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where fi = cos#. 

If polarization has to be accounted for (see JSJ), then the scattering matrix of 
resonant scattering can be written as 



(9) 





where /; , I r , U are the Stokes parameters as defined bv lHamiltonl ( 1947T ) : IChandrasekharl 
( 1950T I. 

For the changes in the surface brightness (see 32J), the difference between the 
isotropic and Rayleigh phase functions is not very important, while it is crucial if 
one needs to evaluate the degree of polarization (see ij5|. Isotropic scattering does not 
produce polarization. Instead it "erases" information on the initial direction and orien- 
tation of the electric vector in the "memory" of the scattered photon. On the contrary, 
Rayleigh scattering changes the polarization state of the radiation field. 
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Fig. 1 Left: Optical depth of the Perseus cluster in the 6.7 keV FeXXV line over the range 
of radii from r to 00. For n ~ 70 kpc the optical depth is of order 1. To a first approximation 
the surface brightness is expected to be only weakly disturbed by resonant scattering for 
projected radii R 2> r\. For smaller projected radii the surface brightness profile is expected 
to be flat. Right: Radial profiles of the He-like iron K a line with (thick solid line) and without 
(dashed line) resonant scattering in the Perseus cluster. Resonant scattering suppresses the line 
intensity in the core and redistributes line photons to larger radii. Pure thermal broadenin g and 
a flat radial abundance profile were assumed in this calculation. Adapted from Chur azov et alj 
J2004D . 



3 Impact on surface brightness profile and abundance 



Since the density of the gas declines faster than ~ r in the outer regions of all 
clusters, the optical depth of the outer regions in any line is small (see Fig. [T] left). 
We therefore do not expect strong modifications of the line surface brightness profile 
at large projected distances from the cluster center (see Fig. [T] right). The central 
part on the contrary is dense and in the region where the optical depth is of order 
unity one can expect that resonant scattering will erase all features in the surface 
brightness, leading to a flattened (s ee Fig. [T] right) surface brightness distribution 
(Gilfan ov. Sunvaev fe Chu razov 1987]). Due to the conservative nature of resonant scat- 
tering, the photons removed from the line of sight going through the cluster center are 
re-distributed to larger projected distances. This effect causes a slight increase of the 
surface brightness at R > r\ (see Fig. [T] right). Since the surface brightness in the 
continuum is not modified, the resonant scattering causes a decrease of the equivalent 
width of the line in the cluster center and a slight increase of the equivalent width in 
the outer regions. This could cause an apparent decrease in the abundance of heavy 
elements in the central part of the cluster if resonant lines are used to measure it 
(|Gilfanov. Sunvaev fe Churazovll 19871 ). 

This effe ct was suggested to play a role in the spurious "abundance holes" in cluster 
centers (e.g. iBohringer et al.l 2001: Sanders et al. 2004), although other effects such as 
the presen ce of multi-t emperature plasma may be more important in these environ- 
ments (e.g. lBuotelkoOOr i. 



Note that optically thin lines are not modified by resonant scattering and therefore 
the observed spectrum might have line ratios different from the expected values for an 
optically thin thermal plasma. This implies that one can use the line ratios (the flux 
in an optically thick line divided by the flux in an optically thin line) rather than the 
line equivalent width (i.e. the ratio of the flux in the optically thick line to the spectral 
intensity of the continuum) . The ideal situation is when two lines from the same ion 
are used, since in this case many issues related to the modeling of the ionization state 
can be avoided. 

A set of simple expressions was d erived by iGilfanov. Sunvaev fe Churazovl (| 19871 ) ; 
ISazonov. Churazov fe Sunvaevl (2002) for model gas density distributions and in the 



limit of t <C 1 for the isotropic and Rayleigh phase functions. These expressions can be 
used for making order of magnitude estimates, while for real systems it is necessary to 
solve a radiative transfer problem to predict the shape of distortions, especially when 
substantial temperature gradients are present, which affect the ionization state and 
therefore the radial distribution of ions. 

Overall the resonant scattering has been searched for in a number of objects, al- 
though clear evidence for the effect was found only in few systems - for the Ne-like iron 
(Fe XVII) line at 15 01 A in NGC4636 l|Xu et al.ll2o"ol iKahn et al.ll2003l :| Werner et al 



20091 : lHavashi et al.l l2009) and for NGC1404, NGC5813 and NGC4472 l|Werner et al 



2009). Nevertheless, the magnitude of the effect (modification of the surface brightness 
distribution in the line) is well within the capabilities of modern X-ray telescopes. In 
fact, one can use the absence of resonant scattering effects to infer the properties of 
the gas velocity field, as discussed in the next section. 



4 Sensitivity to the velocity field 

If one accurately measures the temperature of the intracluster medium (ICM), then 
the line ratios of (optically thick) resonant and (optically thin) non-resonant lines can 
be used as a powerful tool to measure the line widths and therefore to get information 
on the velocity field, which is difficult to obtain by other means. Given that the optical 
depth sensitively depends on the turbulent broadening of resonant lines, especially for 
heavy elements like iron, the comparison of the fluxes of e.g. iron and nickel lines could 
be used to place constraints on the level of turbulence in clusters ( Churazov et al.l 



I2004T ) . For the Perseus cluster, the XMM-Newton spectra of the central region do not 
show a suppression of the resonant 6.7 keV line of iron relative to a much more optically 
thin nickel line when the APEC ( Smith et al.ll200lr ) plasma emission model is used, 



although one should expect the 6.7 keV line of He-like iron to be suppressed in the 
core if the iron line were only thermally broadened. Indeed, without additional line 
broadening the optical depth of the 6.7 keV line is ~ 3 (see Table [I}. Therefore, the 
line in the central ~ l' region must be suppressed by a factor of about 1.7 (see Fig. [2J , 
in contrast with observations. 

Even more robust constraints are expected if optically thic k and optical l y thin 
lines of the same element or even better of the same ion are used. IWerner et al.l ( 20091 ) 



obtained high resolution spectra of the giant elliptical galaxy NGC 4636 using the 
grating spectrometers on the XMM-Newton satellite. A detailed study of the spectra 
proved that the Fe XVII line at 15.01 A is suppressed (Fig. [3]) only in the dense core 
and not in the surrounding regions, while the li ne of Fe XVII at 17.05 A is optically 
thin and is not suppressed. IWerner et al.l 1 20091 ) modeled the radial intensity profiles 
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Fig. 2 Left: 5-9 keV spec trum of the 30" to 2' annulus centered on NGC 1275 and fitted with 
APEC llSmith et al,|[2001J ) models. Right: Influence of turbulence on the strength of the reso- 
nant scattering effect. The curves show the expected ratio of the optically thick and optically 
thin lines calculated for effective Mach numbers of 0, 0.25, 0.5 and 1 (from top to bottom). For 
comparison, crosses show the observed ratio of heavy element abundances obtained ignoring 
parts of the spectrum containing the 6.7 keV and 7.9 keV complexes, respectively. Given that 
the resonant K a line of He-like iron contributes about 50% to the 6.6-6.8 keV complex of lines, 
the measured abundance ra tio is consistent with the curves for Mach number >0.5. Adapted 
from lChurazov et al.l 1120041) . 



of the optically thick line, accounting for the effect of resonant scattering for different 
values of the characteristic turbulent velocity. Comparing the model to the data, it was 
found that the isotropic turbulent velocities on spatial scales smaller than w 1 kpc are 
less than 100 km/s and the turbulent pressure support in the galaxy core is smaller 
than 5 per cent of the thermal pressure at the 90 per cent confidence level. 

While there is still a bit of a controversy in the results, which definitely suffer from 
the limited energy resolution of the present-day CCDs or from the finite size of the 
emitting region for RGS data, it is clear that the magnitude of the effect is within the 
reach of the current generation of telescopes and the situation will improve in future. 

The relation between the resonant scattering effects and the velocity field can 
be used to test even m ore subtle effects, such as e.g. anisotropy of gas motions. 
IZhuravleva et al.l (2010bj) showed that for pure radial motions versus pure tangential 
motions the amplitude of the effect can change by a factor of 1.5 for realistic values of 
the velocity amplitude (see the left panel in Fig. [3J. Clearly, tangential gas motions do 
not change the optical depth in the line (the right panel, Fig. 2} calculated along the 
radial direction. Photons coming from the cluster central region are expected to make 
the largest contribution to the scattered flux. Since they move essentially along the 
radius, tangential gas motions do not affect the probablity of these photons to be scat- 
tered. Therefore, to a first approximation the presence of pure tangential motions does 
not affect much the magnitude of the resonant scattering effect. Radial motions, on 
the contrary, strongly affect the optical depth for the photons coming from the cluster 
center (the right panel, Fig. 0]) and decrease the efficiency of t he resonant scattering. 
The net result of the simulations bv lZhuravleva et al.l (2010bJ) is that the conversion 
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Fig. 3 Left: XMM-Newton RGS spectra extracted from a 0.5' wide region centered on the 
core of NGC 4636. The full line indicates the best optically thin single-temperature plasma 
model fit. The 13.8-15.5 A part of the spectrum, where the strongest Fexvn and Fexvill 
resonance lines are present, was excluded from these spectral fits. Right: Simulated radial 
profiles of the ratio of the 15.01 A line intensities calculated with and without the effects of 
resonant scattering, for isotropic turbulent velocities corre sponding to Mach nu mbers 0.0, 0.25, 
0.5, and 0.75 and a flat abundance profile. Adapted from Werner ct al. (2009). 



of the observable magnitude of the flux suppression to the characteristic amplitude 
of gas motions depends on the anisotropy of the velocity field. In many situations 
the assumption of isotropic stochastic motions is reasonable. Bu t there are cases, e.g . 



spherical shocks/sound waves propagating through the ICM (e.g. lForman et al.ll2007l ) 
wh en strong anisotropy of th e gas motions is expected. 

IZhuravleva et al.l ( 2010bJ) also tested the impact of the spatial scales of motions on 
scattering, showing that suppression of the line flux towards the center of the cluster 
is insensitive to the presence of gas motions on large scales. This is easy to understand 
since small scale motions directly affect the line width and the optical depth of the line. 
If gas motions are present on scales larger than the size of the central region where 
r is of order unity, then these motions can be interpreted as a motion of the whole 
central region. The optical depth and the magnitude of the flux suppression remain 
unaffected. 



5 Polarization due to resonant scattering 



Scattering in strong resonant lines that have a non-zero weight of the Rayleigh phase 
function can lead to linear polarization of scattered radiation (see section I2.2|l . To 
produce net polarization during the Rayleigh scattering of initially unpolarized light 
a quadrupole component has to be present in the radiation field. It appears quite 
naturally in galaxy clusters due to (i) the centrally peaked gas density distribution and 
(ii) differential gas motions (see Fig. [5J. 

In the cluster center, the radiation field is isotropic (for a spherical cluster). There- 
fore, due to the symmetry radiation will be unpolarized. However, at larger distances 
from the center a quadrupole component appears in the radiation field due to the strong 
dominance of the flux coming from the cluster center compared to all other directions. 
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Fig. 4 Left: Impact of the anisotropy of stochastic gas motions on the ratio of the fluxes in 
the Fe XXV line at 6.7 keV in the Perseus cluster calculated with and without allowance for 
scattering. The RMS of the 3D velocity field is 500 km/s in all cases (except for the case with 
no gas motions). Adapted from Zhuravlcva ct al. (2010b). Right: Optical depth of the Perseus 
cluster in the 6.7 keV FeXXV line over the range of radii from r to oo, calculated assuming 
different directions of stochastic gas motions. 



These photons coming from the center give rise to the polarization of the scattered 
radiation (see Fig. [5] the left panel). Obviously, the polarization plane is perpendicular 
to the radius. 

The polarization degree was fi rst calculated bv lSazonov. Churazov fe Sunvaevl ( 20021 ) 
(see also IZhuravleva et al.ll2010l ) for the Perseus and Virgo clusters (see Fig. [6]). For 
the Perseus cluster, the polarization degree was calculated in the He-like iron line at 
6.7 keV, which has an optical depth ~ 3. The polarization is zero at the center of the 
cluster (as expected from the symmetry of the problem) and increases rapidly with 
distance from the center, reaching ~ 10 per cent at ~ 1 Mpc (Fig. [6] the left panel). 
The polarization degree in the Virgo cluster in the most promising lines is a few per 
cent (Fig. [6] the right panel) . 

The assumption of spherical symmetry is a strong one, since clusters are rarely 
spherical. Depending on the characteristics of the substructure in the azimuthal density 
and temperature distributions the polarization degree can change either way: increase 
or decrease even if one neglects non-zero gas velocities. When working with simulated 
clusters we see variations of the polarization degree by as much as factor of two when 
3-dimensional density and temperature distributions are replaced with the azimuthly 
averaged values. 

Gas motions affect both the amplitude and the direction of polarization. If a gas 
lump is moving relative to the whole cluster with high speed, such that the Doppler 
shift of the line exceeds the line width (i.e. the line leaves the resonance), then scat- 
tering will not occur along the direction of motion, while photons coming from the 
perpendicular direction will be scattered, producing net polarization. The sketch of 
one of the possibilities is shown of Fig. [S] (the right panel) . 

IZhuravleva et al.l (120101) . used full 3D models of galaxy clusters taken from cosmo- 
logical simulations (|Dolag et al.l 120081 ; ISpringel. Yoshida. fe White! 120011 ) to calculate 
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Fig. 5 Sketch of the problem. Polarization arises due to the quadrupole component in the 
unpolarized radiation field produced by a centrally peaked gas density distribution (the left 
panel) and/or differential gas motions (the right panel). The line of sight is perpendicular to the 
picture plane. Left panel: scattering occurs at the edge of the cluster (the black dot), the black 
wavy arrows show photons coming from different directions and participating in scattering; 
the red arrow shows the direction of polarization of the scattered radiation (the orientation of 
the electric vector). Right panel: the black dot shows the place where scattering occurs, the 
black straight arrow shows the direction of motion of the gas lump with velocity V relative to 
the whole cluster, the black wavy arrows show photons participating in the scattering while the 
white wavy arrows show photons that will not be scattered due to gas motions (the line leaves 
the resonance in the direction of motion), the red arrow shows the direction of polarization of 
the scattered radiation. 



the polarization degree with an accoun t for gas motions. For one of the most massive 
simulated clusters (g8 in the sample of iDolag et al.l IJ2008T I) the polarization degree in 
the He-like iron line reaches ~ 25 per cent at a distance of ~ 500 kpc from the center 
if no turbulent motions are present (Fig. [7] the left panel). Inclusion of gas motions 
substantially decreases the polarization down to ~ 10 per cent (Fig. [7] the right panel) 
within the same region. They als o made calcu l ation s for a bullet-like cluster ("g72" 
cluster in the simulated sample of Dolag ct al. (2008)), showing that the polarization 
degree is less that 7 per cent within the region where interaction of two sub-clusters 
occurs. 

When polarization in lines is considered, one should take into account the contam- 
ination of the polarized flux in the resonant line by the unpolarized emission in the 
neighboring spectral continuum and nearby lines. The effect is severe. For example, 
for the 6.7 keV line an energy resolution of ~50 eV leads to a factor of 2 drop in the 
degree of polarization. Also, at large distances from the cluster center unpolarized CXB 

emission starts to dominate. 

C urrently new concep t s of X-ray polarim eters are under discussion (see lMuleri et al.l 

120091 : ISoffitta et al.ll2o"o"sl : ICosta et al.ll200sl . e.g.). Galaxy clusters are promising but 
challenging targets for future polarimetric observations. Polarimeters should have good 
energy resolution (better than 100 eV), which is hardly reachable with the nearest fu- 
ture X-ray polarimeters. In addition, arcminute angular resolution is needed to resolve 
nearby clusters. Naturally, a la rge field of view and eff ective area are required to collect 
large numbers of photons (see IZhuravleva et al.ll2010l ). 



6 Impact on the line shape 

The optical depth given by eq. [7] corresponds to the center of the line, i.e. to a photon 
energy equal to the line energy Eq. If a photon is born in the wing of the line, then 
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Fig. 6 Left panel: Expected polarization degree as a function of projected distance from 
the center of the Perseus cluster in the most prominent resonant line of Fe XXV at 6.7 keV 
(calculated assuming constant iron abundance). Right panel: Expected polarization degree 
as a function of projected distance from the center of the Virgo/M87 cluster in the most 
prominent resonant lines. 
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Fig. 7 Polarization degree of the simulated massive cluster g8 l lDolae et al.ll20Qgh in the K a 

line of Fe XXV at 6.7 keV. The polarization degree was evaluated as P = \J Q 2 + U 2 /I. I is 
the total intensity, including scattered and direct emission. The colors in the images denote 
the polarization degree in per cent. The short dashed lines show the orientation of the electric 
vector. Contours (factor of 4 steps in intensity) of the X-ray surface brightness in the chosen 
line are superposed. The size of each picture is 2x2 Mpc. The left panel shows the case of gas 
being at rest. The right panel shows the case of the gas velocities obtained in the simulations. 
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Fig. 8 Simulated spectra of the 6.7 keV line emerging from the core of the Perseus cluster 
(left) and at a projected distance of 800 kpc (right). Lines are broadened only due to thermal 
motions of ions. The calculations were done with (black curves) and without (blue curves) 
accounting for resonant scattering. The line in the cluster outskirts (right panel) is broader 
due to the higher temperature in the outer regions of the Perseus cluster. 



the optical depth will be lower, 



t(E) = r e 






(10) 



If t(E) <C 1, then such photons will leave the cluster without scattering, while photons 
with E ~ Eq will be scattered from the line of sight going through the cluster center. 
It is therefore clear that resonant scattering should affect the shape of optically thick 
lines. The spectral profiles of the He-like iron line at 6.7 keV calculated in the core 
of the Perseus cluster and at a projected distance of ~ 800 kpc are shown in Fig. [H] 
The calculations were done with (black curves) and without (blue curves) allowance 
for scattering. The lines are broadened only due to thermal motions of ions. One can 
notice the suppression of the emission in the line center (left panel) caused by resonant 
scattering: photons with the energy close to the line central energy are scattered from 
the line of sight because of the large optical depth, while photons born in the wings 
escape freely. At the edge of the cluster the optical depth is small and the line profiles 
almost match each other. 



7 Various extensions 



7.1 Cosmological tests 



The effects of resonant scattering are proportional to the optical depth of the cluster 
(see eq. [SJ), which scales linearly with the size of the cluster I and the gas density 
n e : t oc n e l x F(T,Z), where the factor F(T,Z) accounts for the element abundance, 
ionization equilibrium and the line broadening. If another characteristic of the cluster 
(with a different dependence on I and n e ) is known, then both I and n e can be deter- 
mined simultaneously. For instance, one can use the X-ray surface brightness, which is 
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proportional to n e l x W(T, Z), where the W(T, Z) factor accounts for the temperature 
dependent emissivity and the gas metallicity. Once I of a spherically symmetric cluster 
is known, one can immediately perform an angular diameter-redshift test, analogously 
as is done using a combination of the cosmic microwave background distortions (the 
Sunyaev-Zel'dovich effect) and the X-ray surface brightness. Thus a simple cosmolog- 
ical test based solely on X-ray data is possible. This is an attractive possibility since 

such a test can be performed by a single X-ray telescope in a sing le observation. 

Se veral flavors of the same technique have been proposed. iKrolik fc Raymond! 
(|1988l ) and lSarazinl (1989) suggested to look for absorption lines in the spectrum of an 
X-ray emitting background quasar. The equivalent width of the strongest absorption 
lines can be of order several eV. Another po ssibility is to use the distortions of the 
cluster emission itself, u sing the polarization ( SjaaQJlQX^-Ciiuraz pv fc SunvaevH2002r ) or 



shape of emission lines ( Molnar. Birkinshaw. fc MushotzkvM2006r ). or directly the ra- 



dial variations of the equivalent width or the flux ratios o f optically thick and thin 
lines (as discussed by iGilfanov. Sunvaev fc Churazovlll98'7l ). Some variants of these 
methods are within the reach of currently operating X-ray telescopes, but the un- 
certainties in line widths and the complex spectra of the central regions of cool core 
clusters complicate the analysis. Significant progress is expecte d with the l aunch of 
X-ray micro-calorimeters with energy resolution of a few eV fe.g. lMitsuda et al]l2008h . 
They will simultaneously constrain the presence of stochastic gas motions and allow 
for accurate measurements of line shapes (see e.g. Fig[SJ). 



7.2 AGN X-ray echo in lines 

To better understand how massive black holes grow in galactic nuclei, one would like 
to have information on AGN variability on time scales much longer than the several 
tens of years probed by historic light curves. While the behavior on the longest time 
scales (hundreds of millions of years) is determined by the supply of gas (e.g. triggered 
by a merger) from the outer regions of a galaxy to the central kiloparsec, AGN can 
switch on and off on intermediate time scales for a number of reasons, e.g. due to 
instabilities in the accretion disk or due to feedback of the growing massive black hole 
on the surrounding medium. Clearly, such long-term phenomena can only be studied 
by means of indirect observations. An attractive idea is to observe X-ray radiation 
which was emitte d in the past by a galactic nucleus and l ater scattered in our direction 



by ambient gas l|Sazonov. Sunvaev. fc Cra mphorn 2002). Cluster dominant elliptical 



galaxies are the most promising targets for such observations, since their extended gas 
haloes make it possible to probe time scales up to a few million years. Since AGN 
radiation is much harder than the thermal emission of intracluster gas, one possibility 
is to search for AGN "echoes" at energies E 3> kT/(l + z) (T being the gas temperature 
and z the redshift of the object). An additional significant advantage is provided by 
resonant X-ray lines, since the scattered/intrinsic flux ratio is expected to be larger by 
a f actor of 3-10 in a resonance line than in the neighboring continuum. 

ISazonov. Sunvaev. fc Cramphornl ( 2002T ) assessed the level of constraints that could 



be derived from future observations on the past X-ray luminosity of M87 and Cyg A 
(see Fig.[9|, giant elliptical galaxies residing in the central regions of two nearby clusters 
of galaxies. For instance, scattered line radiation should be detectable from the Virgo 
cluster if the X-ray luminosity of M87 was a few 10 erg s _ (i.e. three orders of 
magnitude brighter than presently) until a few 10 years ago. The same method can 
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Fig. 9 Left: Illustration of the delayed X-ray line echo of an AGN outburst in a galaxy 
cluster. Locations with equal travel time from the AGN to the observer lie on a paraboloid. 
In the switch-off case (the AGN was bright for a long time until some instant in the recent 
past), the cluster plasma located in the grey and black regions is scattering the X-ray radiation 
emitted by the AGN. In the case of a short AGN flare, only the grey area is filled by photons 
and contributing to the scattered radiation. A possible path for a photon scattered into the 
observer's line of sight is marked by a dashed line. Right: M87/ Virgo surface brightness of 
scattered AGN radiation (solid lines) in the 6.70 keV line of He-like iron at times (from top to 
bottom) of 5, 50, 100 and 250 thousand years after an AGN switch-off. The dashed line shows 
the intrinsic surface brightness profile of the cluster gas in the 6.70 keV line. The AGN lumi- 
nosity was Lx = 0.01 L^dd before switch-off. Adapted from Sazon ov. Sunvaev. fc C ramphorn 
J2002T) . 



also be applied to groups of galaxies and isolated giant elliptical galaxies, which too are 
large reservoirs of hot gas capable of scattering AGN radiation. Of course, the testable 
time scales are shorter (up to a few times 10 years) than for clusters of galaxies. 



7.3 Warm-Hot Intergalactic Medium (WHIM) and blazars 



Resonant scattering may be important not only in the dense gas inside a cluster, 
but also in a much more tenuous warm intercluster medium in filaments. In the local 
Universe such structures have linear sizes of 5-10 Mpc, temperatures of ~ 10 to ~ 10 
K and density 8 — 10—100 times larger than the mean baryonic density of the Universe 
(Pbar — 3.6 10 g cm ). Simulations suggest that these structures may account for 
a sig nificant or even dominant fraction of the baryons at z — (e.g. ICen &: Ostrikerl 

1 199m . 

The low density of filaments makes direct detection of their thermal emission ex- 
tremely difficult. If there is a bright QSO t hen one can try to se arch for absorption 
lines using high resolution spectroscopy (e.g iNicastro et al.ll2005l ). although it is not 
clea r if existing data have already provided a robust detection of truly diffuse gas 
(e.g Rasmusscn et al. 2007). If one is looking for WHIM in the directions where no 
bright QSO is present, then one can still use a combined signal of numerous back- 
ground sources (and search for an absorption feature) or instead look "between" the 
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Fig. 10 Left: The X-ray spectrum of warm gas with an overdensity of <5 = 30, a temperature 
of 10 6 K, a Thomson optical depth of 10 — 4 and a metallicity of 30% solar. These conditions 
should be typical for a filamentary structure in the warm IGM. The dotted curve shows the 
emission spectrum due to collisional excitation and ionization. For the green curve photoioniza- 
tion by the CXB was taken into account. The black solid curve shows Thomson plus resonantly 
scattered radiation. The upper solid line is the intensity of the CXB from the same region. The 
lower solid line shows 10% of the CXB intensity (the level to which discrete sources may be 
removed). The spectra were convolved with a Gaussian with a FWHM of 2 eV. Right: The 
upper curve shows the total spectrum of a region of the sky containing a filament with Thom- 
son optical depth of 10 — 4 . This spectrum is a sum of the thermal emission of the filament, 
scattered CXB and emission of the background sources, part of which is either absorbed or 
scattered by the filament. Note that resonance scattering does not change the CXB flux and 
only emission lines due to the thermal emission of the filament are visible in this spectrum. 
The lower curve shows the spectrum of the same region if 90% percent of the background is 
resolved into compact sources and removed. Note that the most prominent OVII line at 0.57 
keV becomes visible only in this spectrum. In the hypothetical case that all CXB is resolved 
and removed the residual spectrum should be equal to the sum of the thermal emission of 
the filament a nd scattered CXB emi ssion. These two components are shown in the left panel. 
Adapted from lChurazov et aO (l200lh . 



background sources. As pointed out bv lChurazov et al.l ( 200 if ) the scattering of X-ray 
background photons in He and H-like ions of heavy elements can exceed the "local" 
thermal emission of a filament by a factor of a few or more. Due to the conservative 
nature of resonant scattering, this resonantly scattered radiation can only be identi- 
fied if a significant fraction of the CXB is resolved and removed. While the combined 
spectrum of the resolved sources will contain X-ray absorption features, the residual 
background will contain corresponding emission features with the (on average) same 
intensity. At the relevant densities and temperatures, the lines of He and H-like oxygen 
at 0.57 and 0.65 keV are most promising. These lines (which have a typical width of 
~ 1-2 eV) may contain up to 50% of the total 0.5-1 keV emission of the filament 
(Fig, imp . Note that the detection of WHIM in absorption using a combined signal of 
background objects and the detection of a line in emission using a part of the image 
without bright sources impose different requirements on the telescope characteristics. 
In the former case, the energy resolution is crucial, since the goal is to detect a weak 
absorption line in the combined spectrum of all resolved sources - in practice the best 
results will be achieved if the energy resolution is better than the equivalent width of 
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the absorption line. In the latter case, the same total flux is compared with the unre- 
solved background - the low internal background is crucial in this case, while superb 
energy resolution is less of an issue. In both cases the angular resolution has to be good 
enough to resolve > 50% of the CXB. 

On average up to a few percent of the soft CXB could be resonantly scattered by 
this phase of the IGM and resonantly scattered photons should account for a significant 
fraction of the truly diffuse background at low energies. Close to bright X-ray sources 
such as galaxy clusters or AGN the flux of scattered radiation will be further enhanced. 
From this point of view, off-line blazars are the most promising illuminating sources. 
The scattered emission from AGN may also constrain the duration of the active phase 
of these objects analogously to what is discussed in H7.2I 



8 Conclusions 

Resonant scattering in the brightest X-ray emisson lines leads to a number of character- 
istic observational features, inluding distortions in the surface brightness, modifications 
of the line spectral shape and polarization. The sensitivity of the resonant scattering 
effects to the gas velocity field makes them a promising tool for IGM studies. Resonant 
scattering can also be used for the basic cosmological angular diameter/redshift test, 
as a tracer of a very rarefied IGM and as a diagnostic of past outbursts from AGN. 
While the distortions in line surface brightness profiles are within the reach of modern 
X-ray telescopes, future high energy resolution and polarimetric missions will be able 
to fully exploit resonant scattering as a diagnostic tool. 

Acknowledgements We are grateful to Klaus Dolag, William Forman, Nail Inogamov, Dmitry 
Nagirner, Leonid Vainshtein and Norbert Werner for useful discussions. 



References 

Acton L. W., 1978, ApJ, 225, 1069 

Akimoto P., Furuzawa A., Tawara Y., Yamashita K., 1999, AN, 320, 283 

Akimoto F., Furuzawa A., Tawara Y., Yamashita K., 2000, AdSpR, 25, 603 

Anders E., Grevesse N., 1989, GeCoA, 53, 197 

Bohringer H., et al., 2001, A&A, 365, L181 

Buote D.A. 2000, ApJ, 539, 172 

Cen R., Ostriker J. P., 1999, ApJ, 514, 1 

Chandrasekhar S., 1950, Radiative Transfer, Oxford, Clarendon Press, 1950 

Churazov E., Haehnelt M., Kotov O., Sunyaev R., 2001, MNRAS, 323, 93 

Churazov E., Forman W., Jones C, Sunyaev R., Bohringer H., 2004, MNRAS, 347, 29 

Costa E., et al., 2008, Space Telescopes and Instrumentation 2008: Ultraviolet to Gamma Ray, 

eds M.J.L.Turner,K.A.Flanagan, proc. of the SPIE, v. 7011, 70110F 
Dolag K., Borgani S., Murante G., Springel V., 2009, MNRAS, 399, 497 
Dupke R. A., Arnaud K. A., 2001, ApJ, 548, 141 
Elwert G., 1956, ZA, 41, 67 
Ezawa H., et al., 2001, PASJ, 53, 595 
Forman W., et al., 2007, ApJ, 665, 1057 
Gastaldello F., Molendi S., 2004, ApJ, 600, 670 

Gilfanov M. R., Sunyaev R. A., Churazov E. M., 1987, SvAL, 13, 3 
Hamilton D. R., 1947, ApJ, 106, 457 
Hayashi K., Fukazawa Y., Tozuka M., Nishino S., Matsushita K., Takei Y., Arnaud K. A., 

2009, PASJ, 61, 1185 



18 



Kaastra J. S., Bleeker J. A. M., Mewe R., 1999, NuPhS, 69, 567 
Kahn S. M., ct al., 2003, ASPC, 301, 23 
Krolik J. H., Raymond J. C, 1988, ApJ, 335, L39 
Mathews W. G., Buote D. A., Brighcnti F., 2001, ApJ, 550, L31 

Mitsuda K., Kelley R., McCammon D., Herder J., Ohashi T., 2008, Journal of Low Tempera- 
ture Physics, v. 151, 703 
Molcndi S., Matt G., Antonclli L. A., Fiore F., Fusco-Femiano R., Kaastra J., Maccaronc 

M. C, Perola C., 1998, ApJ, 499, 608 
Molnar S. M., Birkinshaw M., Mushotzky R. F., 2006, ApJ, 643, L73 
Mulcri F., et al., 2009, Journal of Instrumentation, 11, 2 
Nicastro F., ct al., 2005, Nature, 433, 495 
Ostriker J. P., Cen R., 1996, ApJ, 464, 27 
Rasmussen A. P., Kahn S. M., Paerels F., Herder J. W. d., Kaastra J., de Vries C, 2007, ApJ, 

656, 129 
Rugge H. R., McKenzie D. L., 1985, ApJ, 297, 338 
Sakelliou I., et al., 2002, A&A, 391, 903 

Sanders J. S., Fabian A. C., Allen S. W., Schmidt R. W., 2004, MNRAS, 349, 952 
Sanders J. S., Fabian A. C, 2006, MNRAS, 370, 63 
Sarazin C. L., 1989, ApJ, 345, 12 

Sazonov S. Y., Churazov E. M., Sunyaev R. A., 2002, MNRAS, 333, 191 
Sazonov S. Y., Sunyaev R. A., Cramphorn C. K., 2002, A&A, 393, 793 
Shigeyama T., 1998, ApJ, 497, 587 

Smith R. K., Brickhouse N. S., Liedahl D. A., Raymond J. C., 2001, ApJ, 556, L91 
Sofntta P., et al., 2008, Space Telescopes and Instrumentation 2008: Ultraviolet to Gamma 

Ray, eds J.L.Turner, K.A.Flanagan, proc. of the SPIE, v. 7011, 701128 
Springel V., Yoshida N., White S. D. M., 2001, NewA, 6, 79 
Syunyaev R. A., 1982, SvAL, 8, 175 
Werner N., Zhuravleva I., Churazov E., Simionescu A., Allen S. W., Forman W., Jones C, 

Kaastra J. S., 2009, MNRAS, 398, 23 
Xu H., et al., 2002, ApJ, 579, 600 
Zhuravleva I. V., Churazov E. M., Sazonov S. Y., Sunyaev R. A., Forman W., Dolag K., 2010a, 

MNRAS, 403, 129 
Zhuravleva I. V., Churazov E. M., Sazonov S. Y., Sunyaev R. A., 2010b, Astronomy Letters, 

in preparation 



